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Abstract 

Although infinite cylinders are not astrophysical entities, it is possible to 
learn a great deal about the basic qualitative features of generation of gravi- 
tational waves and the behavior of the matter conforming such shells in the 
limits of very small radius. We describe the analytical model using kinetic 
theory for the matter and the junction conditions through the shell to obtain 
its equation of motion. The nature of the static solutions are analyzed, both 
for a single shell as well as for two concentric shells. In this second case, 
for a time dependent external shell, we integrate numerically the equation of 
motion for several values of the constants of the system. Also, a brief des- 
cription in terms of the Komar mass is given to account for the gravitational 
wave energy emitted by the system. 



1 Introduction 

Recently there has been a growing interest in the literature to study the gravita- 
tional dynamics of infinite cylindrical shells of counter-rotating particles in dif- 
ferent geometrical settings[|T||. One of the best studied case is that considered 
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by Apostolatos and Thorne|l2l|. They made a detailed analysis of the system us- 
ing C-energy balance arguments[3| and a sequence of momentarily static and 
radiation-free shells to avoid the complications causes by gravitational radiation 
for a fully relativistic cylindrical shell. However, for a shell outside equilibrium 
the radiation-free requirement is incompatible with the shell equation of motion, 
obtained from the matching conditions of the interior and exterior vacuum solu- 
tions. In this paper we show that there are two equilibrium solutions, correspond- 
ing to constants values of shell radius, the rest mass per unit of proper length (A) 
and the angular momentum per unit of rest mass (7) of the particles composing the 
shell. One of the solution is a stable equilibrium while the other is unstable. Next, 
we study a system of two shells, first in a static configuration and show that the 
equilibrium configuration of the external shell may also be stable or unstable in 
a similar way as for only one shell; then, we integrate (numerically) the equation 
of motion of the external shell for different values of the constants of motion and 
initial values of the coordinate radius starting at rest. This integration is possible, 
without having a complete knowledge of the wave solution outside the external 
shell, because the matching conditions and the fact that the internal solutions sa- 
tisfies the time independent field equations. The innermost shell is chosen in a 
static stable configuration, which allows us to have a smooth spacetime geometry 
on the cylinder's symmetry axis. We obtain shells that perform damped oscilla- 
tions, collapse or are locally expanding, depending on the values of the conserved 
quantities and integration constants. We give a brief description of the energy 
emitted by the system in the case of damped oscillations, in terms of a generalized 
Komar mass. 

The paper is organized as follows. In the next section we describe the geo- 
metry and equation of motion of a cylindrical shell of counter-rotating particles. 
This includes the metric and junction conditions and the description of the sur- 
face stress-energy tensor using the kinetic theory in general Relativity adapted to 
the shell. In section 3 we discuss the stable and unstable equilibrium solutions 
of a single shell. In section 4 the dynamics of two concentric cylindrical shells of 
counter-rotating particles is presented and the matter equation of motion is numer- 
ically integrated to show the difl'erent types of solutions, depending on the values 
of the conserved quantities and integration constants. Finally, a brief summary of 
the main results is given in section 5. 
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2 Geometry and equation of motion of a cylindrical 
shell of counter-rotating particles 



2.1 Metric and junction conditions 

We consider a spacetime M = M" U S U with cylindrical symmetry where 
2 is the history of a hollow cylinder composed of counter-rotating particles of 
rest mass equal to unity; M~ (A/"*") is the vacuum interior (exterior) region of the 
cylinder. In the vacuum interior (M") and exterior (A/"*") of the shell, we introduce 
canonical cylindrical coordinates {t,r,z,4>)- The metric takes the form[|3J. 

dsl = e^^^-^^-{dr^-dtl) + e^'^^dz^ + e-^1'-r^d(P^. (1) 

These coordinates are uniquely determined up to the (non-trivial) change of 
scale 

z — <A — ^lA-y". — >e~^r, t — > e~^t. (2) 
The Einstein field equations in the empty space inside and outside the shell are 

1 

+ = (3) 

r 

y,t = 2r(A,.<A,f, 7,r = Kifr^ + iff]). (4) 

Thus, il/(r,t) plays the role of a gravitational field whose static part is analogue 
of the Newtonian potential. The time depended solutions of (|3]) represent gravita- 
tional waves [4] . Equation ([3]) is the integrability condition of equations (0]). The 
coordinates (z,0,r) and the metric function i// are continuous across the shell S, 
while t and the metric function y are discontinuous. Smoothness of the spacetime 
geometry on the axis r = requires that 

7 = and iff finite at r = 0. (5) 

The junction conditions of and M"*" through 2 require the continuity of the 
metric and to specify the jump of the extrinsic curvature K- compatible with the 
stress-energy tensor on the shell. The induced metric on 2 is given by 

dsl = -dT^ + e^^^dz^ + e-^'f^R^dcf)^. (6) 

Here (Ai(t) = if/+{R(T),t+{T)) = i/f_(i?(T),?_(T)); the evolution of the shell is 
characterized by R{t), which is the radial coordinate r at the shell's location and 
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T the proper time of an observer at rest on Z. The Einstein field equations on the 
shell reduce to[j5i| 

Kt-K;^ = ^n(S,j-\sg,j); (7) 
here 5/; es the surface stress-energy tensor, gij is given by ^ [i = (t,z,0)] and 

Ktj = -(Km;a)efey, = (|-, |, ^) ; h', = 6yN,N'^, (8) 

N is the outward unit vector normal to the shell: 

A^ = A^''— +A^'-; A^''>0. (9) 
or ot 

The tangent vector is 

d . d ^ d 

From the normalization conditions (u - u = - I , N ■ N = l)we get 

^ = X± = + Ve-2(r±-^^)+i?2, (11) 

OT 

Similarly, from u- N = 0,we get 

N''=X-, N'=R. (12) 
The calculation of the tensors K- gives the result (we omit the (+) for simplicity) 

Krr = {7,,-ilr,t)R + Xiy,r-i/r,r)-iRX-RX)e^^-^f\ (13) 



^-iJfAX-iffjR 



(14) 



l\R 

K,z = -e^H>l^,rX + ilfjR]. (15) 
All the quantities are computed at the location of the shell S. 

2.2 The surface stress-energy tensor 

To compute the surface stress-energy tensor we shall use the relativistic kinetic 

theory B adapted to the matter within the shell 2. 

The stress-energy tensor at the event x 6 2 is defined in the usual way by 

JP(x) 

= S''^5{n), (16) 
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here f{x,p) is the distribution function; P{x) is the space of momentum p = 
(p"d/dx") in X, corresponding to particles with proper mass m = I, p ■ p = -1, 
and n is the intrinsic element of volume in this space. Introducing a coordinate n 
normal to S, such that n = defines 2, we have x" = (n,x') and 

S'J= [ p-pJfn, (17) 

JP(X) 



n = - — -dp^dp''', g = det(gij) . 

\Pt\ 

Here P(x) = {p : p e P(x),p ■ n = 0}; / is the restriction to 2 of the distribution 
function /. In absence of collisions, / satisfies the Liouville equation in the one- 
particle phase space M = {(x,p) / xeH; p e P(x)}. 

~ ~ df 

L(f)^-f = 0, (18) 
dT 

where (df/dr) means the derivative of f{x\p^) along a phase orbit. SinceO 

S'^ = j_^p'L(f)n , (19) 

we get from (fTSl) and (fT9l ) the conservation law S !j = 0. Similarly the particle 
four-current density, defined by 

N' = [ fp'n , (20) 

verifies 

W.'. = r L(/)7r . (21) 

A simple way of constructing a solution of the Liouville equation is to consider 
/ as a function of constants of the motion. In our case, due to the cylindrical 
symmetry of the source, the components p^ and p^p of the canonical momentum 
are conserved. The component is the particle angular momentum. Thus, a 
solution, f{pz,P^), of Liouville equation, for particles with p^ = 0, and counter- 
rotating is given by 

f(p„p^) = kd(p,)6(pl-J^) . (22) 
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Here J is the modulus of angular momentum per unit mass. We have the same 
number of particles with = 7 and = -7; is a constant related to the density 
of particles at a given point on S. 

The expression (flTl) can be written in the form: 

S'}= r pipU-^'t'^6{p'-)-^\5{p^-J) + 6{p^ + J)]^dp^dp'^ . (23) 

JPix) IPtI 

To compute 5'^ we have to take into account that p'^ = g^'^p,p. The non-zero 
components of S are: 

S** = ^^^Pg**, (24) 

s"-%-^n. (25) 

The quantities rj and p are the surface energy density and pressure. The shell's 
full stress-tensor is 

S = 7] u^u + p{g + u®u-v®v) . (26) 

Here v = e^'^^e^ is the unit vector in the z direction; E = -p-u. 
From p ■ p = - I \N& have: 

'-^=E^-\^w\ ill) 

where w is the modulus of the particle linear momentum per unit mass. Introduc- 
ing the circumference of the shell IttK = 2nRe^^^, we have 

e^^7 7 

w = — = (28) 



Thus 



Using dini) and ([221) we have 



" = ■ (30) 

N = vu , (31) 
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where v is the (surface) number of particles density. Integration of (|20l ) gives 

y=± ^ 2nvR= — ^A, (32) 

where A = (dm/dz) is a conserved quantity: the shell's total rest mass per unit 
Killing length z[|2|. The shell's total rest mass per unit of proper length (dze'^ = 1) 
is : 

A = - N'uiRe-'^dcf) = Ae'^ . (33) 

From (|29l) and (l32l) we have t] = vE. Summing up, the parameters 'R,A,w,E,p,J 
are invariant by the rescaling Q. On the other hand, Ry A are scale dependent. 
The quantity A and J are conserved quantities during a time evolution of the shell. 

2.3 The shell's equation of motion 

From equations (I13I14I15I ) and (|26l ) the junction conditions (|7]) become: 
The zz component: 

i^^,-i^;„ = 4n(p-r]) . (34) 

The (pep component: 

- ^ ~^ + - fAr^ = 47r(;? + 77) . (35) 

The TT component: 

(r - ^)!n - (r - - (A+ - A-) = 4;r(77 + p) . (36) 

Here A = e^'''^~^\RX - RX)\ iff^n = N{tff) is the normal derivative of the gravitational 
wave field. Equations and (|36l ) are equivalent to: 

^:„-^- =-^J4=- , (37) 

Vl 

X'-'X- = -A A Vl+w2 . (38) 

Using the vacuum field equations ([3]) and dH), and the junction conditions (|37]) and 
(1381) . equation (l36l) becomes: 

dr^ l+w^ i?(l+w2)i R(l+w^) 

Equations (I37H39I ) are the same as (24.a-b-c) of Apostolato and Thorne[|2l. 
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Figure 1 : Rest mass per unit proper length for the equilibrium configurations versus the 
linear momentum per unit rest mass of the particles [equation (l46l)1. For w < 0.8836 the 
equilibrium is stable while for w > 0.8836 unstable . 

3 Static solutions 

3.1 The metric 

For a static configuration we have R = 0; R = 0, iff^t = ^,tt = 0. Then, in the vacuum 
outside the shell the field equations ([3]) and © imply 

(/r+(r) = if/^-Kln(r/R) r>R , (40) 
y+(r) = y + i?ln{rlR) r>R, (41) 

where k, iJ/y. and y are constants. Similarly, in the vacuum inside the shell the field 
equations Q and (H)) plus the boundary conditions ^ give 

=H\ r~ = r<R . (42) 

From (fTTI) we have 

X+ = e^^-y- X- = e^^ . (43) 
Then, from (|38l) and (|43l ) we obtain 

g-r = 1-4aVi+w2 . (44) 
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Using that R = 0,i}/^ = 0, i/r „ = and dH]), the equation of motion ([391) takes the 
form: 

where 

A^g^.«(w) = . ^ . (46) 

In the case of an equilibrium configuration of the shell, we have: 

R = « A = A.quuiw) . (47) 

In figure ([B we plot Agg„,7(w); the maximum value Kequii,max = 0.1588 is reached 
for w = 0.8836. Given a constant A < 0.1588 we find from (|46|) two values of 
w that correspond to equilibrium configurations; or, we may use this equation to 
obtain a value of A from w. Recall that w = J/'R; i.e., for fix J, K decreases while 
w increases. For a fix value of w (or A), K is proportional to J. Thus, we may 
have % approaching zero by just taking J very small. 

From the junction condition (|37]) and equations (|40H44I) . (|46l) and (|47]) . we 
obtain 

K = 2w^ ; (48) 
y = ln{\+lw^f . (49) 

Thus, all the scale invariant quantities that define the static configuration are de- 
termined; the scale dependent constant if/j. may be computed from the relations 
g-fe = A//i = "RIR, if a value of /I or i? is given (fixing the coordinate system). In 
the next subsection we analyze the stability of the equilibrium configurations. 



3.2 Stability analysis 



To determine the nature of the static configurations described in the previous sub- 
section we perform a virtual displacement 5% from its equilibrium value. Then, 
we compute from (|45|) the virtual change in the acceleration "R, around = in 
the form 



6K = ^ 

&R 



5% 



(50) 



eq 
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Figure 2: Acceleration vs proper radius of the shell for a momentarily static and radiation 
free shell[2| given by ^5^. The shell with radius Ki - 0.00044 is unstable while = 
0.0024 corresponds to a static stable shell. 

The shell is in a static stable configuration if {dKldK)\ii^^ < 0; similarly the equi- 
librium configuration is unstable if {dKldK)\ii^ > 0. From (|45l) . we have 



From figure ([B we immediately conclude that the stable static configurations co- 
rrespond to w < 0.8836; i.e., for values oi'R> 7/0.8836 (this includes the New- 
tonian limit); while for w > 0.8836, K < 7/0.8836, the equilibrium configurations 
are unstable. In figure Q we plot "/^ as a function of for 7 = 0.001; A = 0.1. 
Notice that at Ki = 0.00044, {&R/&R)\'r, > it is sensible larger than | dK/&R 
for 'R2 = 0.0024; "Ri and 'R2 are the unstable and stable radius, respectively, of the 
shell. 




(51) 
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4 Dynamics of two cylindrical shells of counter-rotating 
particles 

4.1 The model 

We consider two concentric hollow cylinders composed of counter-rotating parti- 
cles which define singular hypersurfaces £ and If of coordinate radius R and R' , 
respectively {R' > R), such that the spacetime M = M*^ U 2 U M~ U 2' U M+. In 
the vacuum interior (M*^), intermediate region (M"), and exterior (M"*"), we use 
canonical cylindrical coordinates {t,r,z,(p) as before (see ([T])). The hypersurface 2 
is considered in a static, stable equilibrium (i?=const). Therefore the region 
and M~ are related through 2 in the same way as in section 3. In general, the 
radius of the singular hypersurface 2' will be time dependent. Thus, the metric 
coefficients in the different regions are given by: 

a) 0<r<R: 



lAo = ^Ax = constant ; yo = • 

R<r< R'(t) : 

il/-{r) = -Kln{rlR) + il/Y, ; 

y~{r) = I? ln{r/R) + y 

K = 2w^; y = lnil+2wY- 

R'(t) < r : 

ip*{r,t^) ; y^{r,t^) (wave solutions of the field equations) 



(52) 



(53) 
(54) 



The matching conditions at rj' = R'{t) require 



lfj\R'{T),f{T)) = riR'iT)) = ifJviT) , 

tff^'(r) = -Kln(R'{T)/R) + il/^ . 



(55) 



The shells' surface stress-energy tensors are (see subsection 2.2) 



S = 
S' = 



rj u®u-\- p{g-\-u®u-v®v) ; 

T] U®U+p{g+U®U-V®v), 



(56) 
(57) 



with obvious meanings for the symbols. 
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Figure 3 : Rest mass per unit proper length of equilibrium configurations of the external 
shell versus linear momentum per unit rest mass of the particles requation (l6n )1. The 
stable configurations correspond to si < s < S2. For s > S2 the equilibrium is unstable; 
si -0.6957,52 = 1.5346. 



4.2 Static solutions 

A static solution corresponds to R' = const. The metric coefficients in M"*" are 



(/r+(r) = ln(r/R') + if/^' 
y+(r) = -K^ ln{rlR') + y" 



(58) 
(59) 



where k' , tff-z' ,y"are constants to be determined. Equations (1381) and (l39l) {R' = 
= i?') become 

(g-r^ - e-^')\R' = -4A'E' , (60) 
^, _ Vl + s\s^ - lw\\ + ^2)] ^ % 



Where 



(1 +2^2)2(1 +2vt;2) 



E' = VrT7, s = ^, y^(R') = y" , y-(R') = y + ln(R'/Rf = y' , (62) 



■(^) ^K,urM^^^R^J') ■ (61) 
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Figure 4: Acceleration vs proper radius of a momentarily static and radiation free exter- 
nal shell[2|. The plot correspond to equation (l64l ). The shell with radius K'l = 0.003744 
is unstable while 'R'2 = 0.009493 corresponds to a static stable shell. 



J' is the angular momentum per unit mass of the particles constituting If ,'R' = 
R'e-'f'^' . Notice that (EB reduces to (1461) for w = 0. From dSSl) we have 



(63) 



In figure ([3]) we show A'^^^^.^{s,'w,R,J') as a function of s for w = 0.4038, R = I 
(the values corresponding to the stable static shell of section 3) and J' = 0.01. 
Considering R = 0, = 0, = -(2w^ /'R')e~'^ , a similar calculation as the 

one done for the interior shell (S), shows that (l39l) becomes 

.. /7?\i+2,.2 (1+2^2)2 

Thus, the static stable configuration of the exterior shell (2') correspond to those 
values of s, such that (5A' Jds) > 0. For w = 0.4038,!^ = 0.002476, J' = 0.01 
we have a stable configuration for 0.6957 < s < 1.5346 [see figures (I3]and|4l) ; s = 
1.05335]. 

To complete the static solution we need the expressions for k' and y" . A 
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Figure 5: Sequence R'{t) of a damped oscillating solution for A' = 11.6070, /' = 
0.01, R'. = 3. 



lengthy but straight forward calculation gives: 

~ 1+8^2(1+^2) + 8^2(1+^2)2 ' ( ) 

4vf4 

/ \TTi^ 1 + 8w2(l + s^)[s\l + 2^2) + 2^2] 
^ ^ = . . o.... . . • (66) 



n') (1+2^2)2(1+2^2)2 

4.3 Damped oscillating solutions 

Let us now consider a simple time dependent model by assuming that the external 
shell is time dependent: r^' = R'ir). We can integrate the equation of motion 
(l39l) for 'R'ir) without knowing explicitly the metric for r > R'{t). Outside the 
external shell we have a wave solution of the field equations Q |4]) that satisfies 
the boundary conditions 

<A^(i?'(T),f+(T)) = <AE'(T), (67) 

r-'(R'(T),t^(T)) = r"(T), (68) 

where i/'s'(t) is given by (1551 ): y"(T) may be obtained from the matching condi- 
tion (1381) at the hypersurface Z'. We specify a particular system choosing an inter- 
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Figure 6: The first oscillations of a siiell witii R'. = 5, A' = 11 .6070, /' = 0.01. It shows 
the asymmetry of the oscillation when it is quite apart from equilibrium. 



nal static stable shell, values of the constants A' and J' and the initial conditions 
R'(0) = R'., R'iO) = 0. From ([55]) we see that the condition (Kir) /K) > 1 is satisfied 
if R < R'{t), V T. In figure© we show a sequence R'(t) of a damped oscillating 
solution, obtained through a numerical integration of (|39l) for A' = II .6070, J' = 
0.01, R' = 3. The internal solution corresponds to the stable static solution of 
section 3{A = 0.l,'R = 0.002476, J = 0.001, = 1). We see that the shell oscil- 
lates until it reaches the equilibrium radius R'^^ = 2.7547. This radius corresponds 
to the stable configuration with A' = A'^^^^.^{s{Tf),w,R,J'), where r/ > 3 x 10^, 
s(Tf) = 1.05335 < 1.53463. The same equilibrium configuration is reached for 
different values of R'^ in a neighborhood of R^^ (we tried R'^ = 5, and 2). 



4.3.1 Generalized Komar energy 

All the knowledge that we have of the external solution is reduced to as the 
boundary of A/"*"; i.e, the matter of the external shell. We may calculate a (time 
dependent) Komar massfT] per unit of proper length in M"*" = M"*" U 2', using the 
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definition 

1 



Mg = j_^e\S'b"--S'5lY{n)dI.a , (69) 

where = ; £ : = const ; 6{n) = -^== (5(r - i?(T)) ; dLa = yp^drd(pdz ; 
V~i^ = y[^re^^^~^^^\ ^ is a constant. Then (|69l ) becomes 

MG = fcJ J 7?'(T)(5;'--5'je^^'-^^'^6/z# . (70) 
Using 5 = 77'(X+)2 + p'(/z+)^'; {h^f = {g+f + (X+)2, we have 



choosing = -1, (|70l) becomes: 

Mg = A'e^'e-^^' (^=)(^ +i?'2^2(r'-^,,)J (72) 

to obtain y' we have to use the matching condition 

X*-X- = -4 A' (73) 

In figure® we show R'(t) and 7?' for t < 0.5, i?^ = 5, and A' and 7' equal to the 
values corresponding to those of figure ©. It is apparent that in the first instances 
of the motion R' grows faster than it decreases. It may be inferred, according to 
order of magnitude estimates using the reduced quadrupole moment of the shell 
per unit of proper length[8], that the gravitational luminosity is larger after the 
shell bounces when it reaches its minimum radius |l9l . It is clear from (|72l) that 
most of the energy radiated away, during the radial oscillation of the external shell, 
comes from the term proportional to i?'^;i.e, from its kinetic energy as it would be 
expected. 



4.4 Collapsing and initially expanding solutions 

An interesting set of solutions of equation (|39l ) corresponds to those with Rj close 
to an unstable equilibrium configuration {R[,gq)- For instance, if R'^ < R[,gq, the 
solution of (1391) is a collapsing shell that hits the internal shell in a time of order 
Tc ~ 10"^ [see figure (|7])]. On the other hand fori?^ > R'^^gq the solution is initially 
expanding, to settle in a damped oscillating solution with a larger stable radius. An 
example is shown in figure ([8]); the unstable equilibrium radius is R[^gq = 1.3660. 
The shell settles in a stable equilibrium radius R'. « 2.3. 
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0.005 0.01 0.015 0.02 1 c 
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Figure 7: A collapsing solution corresponding to 7?^ = 1.3 and R'^^^ = 1.3660; Tc = 0.0222. 
The constant of the model are w = 0.4038, R^l,J^ 0.001, J' = 0.01, A' = 14.59. 




0.1 2 0.3 0.4 0.5 



Figure 8: An initially expanding solution. The parameters are the same as in figured?]), 
R'. - 1.4. 
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5 Final comments 



We have been interested in the relativistic dynamics of cylindrical shells of counter 
rotating particles. This system have been considered before, in some detail, by 
Apostolatos and ThorneO. Our main contribution to the knowledge of the model 
is the following: 

• We have described the matter composing the shells using the relativistic 
kinetic theory adapted to a singular hypersurface. This approach has the ad- 
vantage that the surface energy momentum tensor satisfies the conservation 
laws if the distribution function is a solution of the Liouville's equation. 

• We have analyzed the stability of the equilibrium solutions. We found 
that the increasing branch of the plot of equation (|46l) represents the sta- 
ble equilibrium, while the decreasing part of the plot are the unstable equi- 
librium configurations. This figure was first introduced by Apostolato and 
Thome|l2]| for another purpose. 

• The study of the equilibrium solutions was extended to two concentric shells. 
The results are similar as for only one shell [Compare figures ([HO with ([3l 

m. 

• Armed with the previous results, we analyzed the damped oscillating solu- 
tions of the external shell in the two shells model. These solutions converge 
to a stable static configuration, depending on the values of the constant of 
the motion and initial conditions. 

• The most interesting result was to show that the solutions that start with a 
radius smaller than the corresponding radius of a unstable static solution, 
collapses to hit the internal shell, whose proper radius may be quite small, 
depending on the angular momentum per unit mass of the particles of the 
internal shell. The spacetime that results is free of singularities as long as 
this angular momentum is not null. 

Finally, we think that more information can be obtained from the equation of 
motion of a time dependent shell, about the basic qualitative features of generation 
of gravitational waves and the behavior of the matter conforming such shells, in 
the limit of very small radius. An interesting problem would be to relate the 
motion of the shell with the C-energy[|2l[3l| carried away by gravitational waves, 
while in the interior of the shell we have a superposition of ingoing and outgoing 
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cylindrical gravitational waves that satisfies the boundary conditions ([5]). This and 
related topics may be discussed elsewhere. 
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